A boundary integral method applied to the 3D water coning problem
نویسندگان
چکیده
Often in oil reservoirs a layer of water lies under the layer of oil. The suction pressure due to a distribution of oil wells will cause the oil-water interface to rise up towards the wells. A three-dimensional boundary integral formulation is presented for calculating the steady interface shape when the oil wells are represented by point sinks. Sophisticated integration techniques are implemented in an effort to obtain accurate results. In particular, the efficiency of various integration methods are compared for this problem, including QUADPACK routines, adaptive methods based on the IMT rule, the Kronrod rule, the method of degenerate quadrilaterals, and the Gauss-Rational rule for infinite integrals. Numerical results for various general multi-sink distributions are discussed, as are some further results for the axisymmetric single well problem.
منابع مشابه
Simulation of Water Coning in Oil Reservoirs Using a Corrected IMPES Method
Implicit pressure-explicit saturation method (IMPES) is widely used in oil reservoir simulation to study the multiphase flow in porous media. This method has no complexity compared to the fully implicit method, although both of them are based on the finite difference technique. Water coning is one the most important phenomenon that affects the oil production from oil reservoirs having a water d...
متن کاملBoundary temperature reconstruction in an inverse heat conduction problem using boundary integral equation method
In this paper, we consider an inverse boundary value problem for two-dimensional heat equation in an annular domain. This problem consists of determining the temperature on the interior boundary curve from the Cauchy data (boundary temperature and heat flux) on the exterior boundary curve. To this end, the boundary integral equation method is used. Since the resulting system of linea...
متن کاملThree dimensional static and dynamic analysis of thick plates by the meshless local Petrov-Galerkin (MLPG) method under different loading conditions
In this paper, three dimensional (3D) static and dynamic analysis of thick plates based on the Meshless Local Petrov-Galerkin (MLPG) is presented. Using the kinematics of a three-dimensional continuum, the local weak form of the equilibrium equations is derived. A weak formulation for the set of governing equations is transformed into local integral equations on local sub-domains by using a uni...
متن کاملA Boundary Meshless Method for Neumann Problem
Boundary integral equations (BIE) are reformulations of boundary value problems for partial differential equations. There is a plethora of research on numerical methods for all types of these equations such as solving by discretization which includes numerical integration. In this paper, the Neumann problem is reformulated to a BIE, and then moving least squares as a meshless method is describe...
متن کاملNumerical Solution of a Free Boundary Problem from Heat Transfer by the Second Kind Chebyshev Wavelets
In this paper we reduce a free boundary problem from heat transfer to a weakly Singular Volterra integral equation of the first kind. Since the first kind integral equation is ill posed, and an appropriate method for such ill posed problems is based on wavelets, then we apply the Chebyshev wavelets to solve the integral equation. Numerical implementation of the method is illustrated by two ben...
متن کامل